MA 4 
(262/1959) UDC 621-52:517.512.2 


ACTA 
POLYTECHNICA 


MATHEMATICS AND COMPUTING MACHINERY SERIES No. 3 


JENS G. BALCHEN and ARNE G. BERRE 


A Method for Evaluating the Accuracy in the 
Time Domain Associated with Approximation 
in the Frequency Domain 


Norwegian Contribution No. 3 
Trondheim 1959 


ACTA POLYTECHNICA SCANDINAVICA 


. a Scandinavian contribution to international engineering sciences 


Published under the auspices of the Scandinavian Council for Applied 
Research 
in Denmark by the Danish Academy of Technical Sciences 
in Finland by the Finnish Academy of Technical Sciences, the Swe- 
dish Academy of Engineering Sciences in Finland, and the State 
Institute for Technical Research 


in Norway by the Norwegian Academy of Technical Science and the 
Royal Norwegian Council for Scientific and Industrial Research 


in Sweden by the Royal Swedish Academy of Engineering Sciences, 
the Swedish Natural Science Research Council, and the Swedish Tech- 
nical Research Council 


Acta Polytechnica Scandinavica is divided into the following sub-series: 


Chemistry including Metallurgy Series, Ch 

Civil Engineering and Building Construction Series, Ci 
Electrical Engineering Series, El 

Mathematics and Computing Machinery Series, Ma 
Mechanical Engineering Series, Me 

Physics including Nucleonics Series, Ph 


For subcription to the complete series or to one or more of the sub-series and 
for purchase of single copies, please write to 


ACTA POLYTECHNICA SCANDINAVICA PUBLISHING OFFICE 
Box 5073 


Stockholm 5 
Sweden ' Phone 61 47 53 


This issue is published by 
NORGES TEKNISKE VITENSKAPSAKADEMI 
Trondheim, Norway 


and is designated as NTVA Publication series 1, No. 5, 1959 
Editor: Knut Alming 


| 

} 

- 


MA 4 
(262/1959) UDC 621-52:517.512.2 


ACTA 
POLYTECHNICA 


MATHEMATICS AND COMPUTING MACHINERY SERIES No. 3 


JENS G. BALCHEN and ARNE G. BERRE 


A Method for Evaluating the Accuracy in the 


Time Domain Associated with Approximation 
in the Frequency Domain 


(Received March 11, 1959) 


Norwegian Contribution No. 3 
Trondheim 1959 


} 
| = 


‘ 
- 


A METHOD FOR EVALUATING THE ACCURACY IN THE TIME DOMAIN 
ASSOCIATED WITH APPROXIMATION IN THE FREQUENCY DOMAIN 


by 


Jens G.Balchen and Arne G. Berre*™) 


SUMMARY. 


A method is outlined which makes it possible by 
simple means to estimate the error in the time re- 
sponse arising from a certain inaccuracy between an 
exact and an approximate system in the frequency 
domain. A new error criterion which makes use of 
eee amplitude and frequency plots is intro- 

uced, 


1. INTRODUCTION 


One common way to describe the behaviour of linear systems 
in control and electronic engineering, is by means of frequency 
response characteristics. These characteristics give the ampli- 
tude and phase relationship between a sinusoidal input signal and 
the resulting (also sinusoidal) output signal. The information 
one can get from an amplitude—phase-—frequency plot will in general 
be sufficient to characterize the system performance. The 
frequency response function could be derived either by an analyti- 
cal procedure or by graphical or experimental methods. 

Even though the information contained in the frequency plot 
gives the full idea about the dynamic behaviour, it is often 
desirable to be able to indicate the response of a known transient. 
In the literature there are several methods for determining the 
transient response when the amplitude—phase-—frequency relationship 
is known [1,2]. Unfortunately most of these methods are compli- 
cated and time consuming and restricted to certain classes of in- 
put functions. The accuracy associated with these procedures also 
is difficult to state. A logical procedure, that is often used in 
practice, is to approximate the actual amplitude—phase—frequency 
plot with a known and easily synthesized characteristic of such 
form that the transient could readily be evaluated. 


2, THEORY 


In this paper it is intended to derive the accuracy neces— 
sary for the approximation in the frequency domain to achieve a 
certain maximum error specified in the time domain. 


x) Automatic Control Laboratory, Norwegian Institute of 
Technology, Trondheim, Norway. 
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Symbols used 


x(t), X(s) -— input signal, Laplace transform 
y(t), Y(s) - output response, Laplace transform 
H(s) - transfer function 


symbols used without indicies refers to actual system. 
Index 1 refers to the approximate system. 


X(j), Y(j@) etc. is used for X(s), Y¥(s), for s=j® etc. 
amplitude ratio 

(jw) — phase shift 

a(t) = y(t) - y(t) - error in time domain 


Fig. 1 


Fig. 1 shows an example of the three quantities from the 
last expression above. It is obvious that the approximation is 
satisfactory when 2(t) has a small magnitude. A reasonable 
measure of the error of the approximation is thus given by the 
following expression 
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A small average walue of |s(t)| during the settling periode conse- 
quently means a good approximation. The integrals in Eq. (1) 
however, do not easily lend themselves to calculation from the 
frequency response characteristics only, but it is possible to ob— 
tain a good approximation to integrals of the form 


S\s(+) lat 
° 

by using the expression 
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and choosing @, and a, so as to maximise the integral. 
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It is easily shown that 


= |¥(jo) - y, (jo) | (3) 


This is demonstrated in more detail in ref. [3]. 
A similar expression can be derived for the denominator in Eq.(1), 


Hence 
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If the frequency characteristics for the system and the input sig- 
nal are known, the expression in Eq. (4) can be evaluated. When a 


polar plot is used, as shown in fig.2, ly(5~) | pax and 


ly( jw) ¥,(3@) | are readily determined. The polar plot, 
however, is not very practical to use when synthesizing the 
approximate system. The Bode diagram. using logarithmic scale 
for amplitude ratio and frequency and linear scale for phase 
shift. is more flexible for this purpose. Because the difference 
in Eq. (8) can not be derived directly from separate Bode 
diagrams of Y(jw) and Y,(j®), Bq. (8) is rewritten 
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The factor 


H 
lp| = |1 - | od (% ~ (6) 


can now be determined graphically, when 
= and AY = 9,- 


are given. From the Bode diagram ln| and lq, | are known in db 
(decibel) and consequently Aln| is immediately given in db. 
Fig. 8 shows |p| (db) as a function of Aln| (ab) and Ad (degrees) 
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Fig.3 


This plot is quite similar to the Nichols chart used for servo— 
mechanisms calculations. Instead of using Bq. (1), the accuracy 
of the approximate system will be formulated 


ly(3~) | 


Eq. (7) can be interpreted in another way, i.e. that the maximum 
error |y(jw) - (30) in the frequency domain is a fraction 


In Eq. (7) ly(5) | nex and € are supposed to be known. 
lx(3~) | nex is known from the amplitude ratio characteristic and 


€ is the accuracy that is required in the time domain. 
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If 1% is specified for the accuracy in the time domain, Eq.(1) 
will give € = 0,01 and accordingly € = ~ 40 db as defined in 
Eq. (7). From Bq. (6) and Eq. (7) it follows 


lel. 


max(4b) = + €(ad) = ab) 


(8) 


From this equation it is possible to determine corresponding 
values for Alx| and AP for all frequencies. As an example the 
following transfer function is chosen 


T,8 
(1+7,8) (1+T os) 

1 
Input signal: X(s) == (step function) 
For Ty = 5 and % * 1 the response y(t) is shown in fig.4. 


Fig. 4 
It is easily shown that 


= 7, = 5 214 ab 


The frequency characteristics for the system chosen is 
plotted in fig. 5. With an € = -40 db it follows that 


(ad) (2d) + E(ab) |y|(av) = -26 - |y| (av) 


The figures for ly( jw) | are found in the amplitude characteristic 
in fig. 5, and by means of fig. 8 the corresponding values of 
and are determined. In fig. 6 the allowable errors Al 
and AP are shown as functions of frequency. It appears that it 
is essential to have small deviations for frequencies where 
x(jo)| is large in order to maintain accuracy, in this case for 
low frequencies. 
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Further it is discovered that the accuracy is independent 
of whether the error in phase shift is positive or negative. 


-4 
a 
“12 
A 
agl2 
4 
-20° 
-40° \ 


ica 
mma 
-5 
40 A 
Fig. 5 
ee 
BS 
be 
Ba 
we 
ME 
O1 02 2 #46 10 
Fig. 6 
6 


As to error in amplitude ratio, negative db-values of Ala| are 

less significant than positive. The preceding discussion also 

clarifies that the input signal plays an important part in the 

accuracy resulting from the approximation in the frequency domain. 
If the response y(t) is such that 


@ 

Jy(t) at 


there is a need for defining the accuracy in another way than is 


done in Eq. (1) and Eq. (7). Such a case would arise when the in- 
put is a step, X(s)=<, and 


H(s) ~ const # 0 


It seems then logical to extend the definition given in Eq. (1) 
to 


a(+)| at 


(9) 
ly(t) - (t) lat 


where y,(t) is the straight line which y(t) approaches when t ~ ™, 
fig.7. 
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Fig. 7 


For the function shown in fig.4. y, (t) = O and accordingly Eq. (9) 
leads to Eq. (1). The modified form of Bq. (7) is 


ly( jw) Yo( 3%) | 


where Y, (jw) is the Laplace transform of yq~ (t+) with s = jo. 
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The expression in the denominator in Eq.(10) can be evaluated in 
the same manner as described in Eq.(5). Doing this, however, 
requires some additional work and it seems reasonable instead to 
make an estimation based on the denominator in Eq. (9). If the 
transient response is known approximately, the area involved 
should be rather easily determined with a sufficient accuracy. 

Another approach is to choose the input function x(t) so 
that y(t) = 0. 

If x(s) =+ 


and H(s) ~ const # O 


X(s) can be modified through the substitution 


choosing the parameter 7 such that 


foe) 
Jy(t) at 
° 
is minimized without the loss of characteristic features associ- 
ated with the original function. 
8. APPLICATION 


In order to demonstrate some of the properties of the 
method outlined in this paper, the system shown in fig. 8. is 


chosen. 
r(t) c(t) 
Ris) (1+ 10s)(1+2s)* C(s) 


r(t) =1 
R(s)s 
s 


Fig. 8 


An exact calculation of c(t) for a given r(t) is very 
complicated. By means of the Nichols chart 4 (jw) is evaluated. 
To approximate c(t), 4 (j) is approximated by functions which 
more easily leads to the time response. In fig. 9. is shown how 
R (j@) is approximated with a second order-system in combination 
with a transportation lag. The parameters of the approximate 
System are chosen so that 


H, (30) 


for jw~~ 0 
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The effect is that 
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The 13) | is evaluated by means of fig. 8. and the 
graphical approach gives a final result 


ax % -18db = 0,22 


This value can be compared with an area equivalent to that shaded 
in fig. 7. With a rise time of about 1,83 sec, an overshoot of 
about 40 o/o, a transportation lag of 2 sec. and yy, (t) = 0,7, 
the approximate area is (0,7 + 0,7 ° 0,42) (2+%° 1,8) % 2,7 


Hence the accuracy 0.22 + 1400 x 8 o/o 


In fig. 10 the exact and approximated responses sare shown. The 
solutions were obtained from an analog computer setup. 
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